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' We have used WFPC2 on the Hubble Space Telescope to obtain F450W and F814W images of four cD 

Cn ! galaxies (NGC 541 in Abell 194, NGC 2832 in AbeU 779, NGC 4839 in AbeU 1656 and NGC 7768 in 

O ■ Abell 2666) in the range 5400 < cz < 8100 km s'^. For NGC 541, the HST data are supplemented by 

ground-based B and / images obtained with the FORSl on the VLT. We present surface brightness 
and color profiles for each of the four galaxies, confirming their classification as cD galaxies. Isophotal 
analyses reveal the presence of subarcsecond-scale dust disks in the nuclei of NGC 541 and NGC 7768. 
O ' Despite the extreme nature of these galaxies in terms of spatial extent and luminosity, our analysis 

of their globular cluster systems reveals no anomalies in terms of specific frequencies, metallicity gra- 
C/3 i dients, average metallicities, or the metallicity offset between the globulars and the host galaxy. We 

' show that the latter offset appears roughly constant at A[Fe/H] ^ 0.8 dex for early-type galaxies 

spanning a luminosity range of roughly four orders of magnitude. We combine the globular cluster 
. . metallicity distributions with an empirical technique described in a series of earlier papers to investi- 

/\^ ' gate the form of the protogalactic mass spectrum in these cD galaxies. We find that the observed GC 

metallicity distributions are consistent with those expected if cD galaxies form through the cannibal- 
ism of numerous galaxies and protogalactic fragments which formed their stars and globular clusters 
before capture and disruption. However, the properties of their GC systems suggest that dynamical 
friction is not the primary mechanism by which these galaxies are assembled. We argue that cDs 
instead form rapidly, via hierarchical merging, prior to cluster virialization. 

Subject headings: galaxies: elliptical and cD — galaxies: individual (NGC 541, NGC 2832, NGC 4839, 
NGC 7768) — galaxies: star clusters — galaxies: structure 



1. INTRODUCTION 

In spatial extent and luminosity, cD galaxies are in- 
termediate between normal galaxies and galaxy clusters. 
Indeed, the cD classification itself was introduced in the 
work of Matthews, Morgan & Schmidt (1964) to denote 
the very large D galaxies that they found in some clus- 

^ Based on observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS 5-26555 

^ Based in part on observations obtained at the European South- 
ern Observatory, for the VLT program 68.D-0130(A). 

^ Claudio Anguita Fellow 



ters. (The 'c' prefix was borrowed from the notation for 
supergiant stars in stellar spectroscopy.) Although the 
extraordinary sizes of cD galaxies was readily apparent 
to Matthews et al. (1964), the outstanding nature of 
their distinctive envelopes was fully appreciated only af- 
ter the deep photographic surface photometry of Oemler 
(1976) and Schombert (1986; 1988). 

In this paper, we adopt the definition of a cD as given 
in Schombert (1986) and Tonry (1987). That is, cDs are 
elliptical galaxies with shallow surface brightness profiles 
dlog^y/dlogr « — 2 (1) 
at fly ~ 24 mag arcsec"^. These galaxies exhibit a 
characteristic 'break' over an r^/** law, beginning at 
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rb ~ 10 — 100 kpc, and arc much brighter than typical 
elliptical galaxies, with luminosities ^ lOL* (Sandagc & 
Hardy 1973; Schombert 1986). cD galaxies are always 
found in dense regions, and in virtually all cases, they 
are located near the spatial and kinematical center of 
their host cluster, or subcluster. These facts would seem 
to suggest that the formation of cD galaxies is unique 
to the cluster environment and is closely linked to its 
dynamical history. 

There are, in fact, reasons to believe that the envelopes 
themselves are distinct entities from the galaxies them- 
selves. First, cD envelope luminosity is weakly correlated 
with some properties of the host cluster, most notably 
with cluster richness and X-ray luminosity (Schombert 
1988). Second, the envelopes have surface brightness pro- 
files with power-law slopes that are similar to those mea- 
sured from the surface density profiles of the surrounding 
cluster galaxies. Finally, both the position angle and cl- 
lipticity of cD galaxy isophotes commonly show discon- 
tinuities at Vb, where the envelope begins to dominate 
the surface brightness profile (Schombert 1988; Porter, 
Schneider & Hoessel 1991). 

Models for the formation of cD galaxies have tended 
to focus on dynamical effects that, in dense environ- 
ments, might produce extended low-surface-brightness 
envelopes. A tidal stripping origin was first proposed 
by Gallagher & Ostriker (1972), and examined in detail 
by Richstone (1976). In this scenario, close encounters 
between galaxies liberate material which surrounds the 
central galaxy; if the velocity of the encounter is suffi- 
ciently low, the interaction may result in a merger. The 
most massive galaxies would be preferentially depleted 
as they are most strongly affected by dynamical friction. 
If many such mergers occur, then the central galaxy may 
grow primarily through this process of "cannibalism". 
The rates for galactic cannibalism in rich clusters were 
computed first with approximate theories that used av- 
erage interaction rates without considering the actual 
galaxy trajectories (Ostriker & Tremaine 1975; White 
1976), and later with Monte Carlo techniques (Hausman 
& Ostriker 1978). These calculations predicted rather 
high cannibalism rates in rich clusters, leading to a lu- 
minosity growth of the first ranked of ^ lOL* over the 
lifetime of the host cluster. The finding that ~ 50% of 
brightest cluster galaxies (BCGs) contain multiple nu- 
clei {e.g., Hoessel & Schneider 1985) has sometimes been 
interpreted as supporting evidence for this model. 

Richstone & Malumuth (1983) carried out more de- 
tailed simulations of the growth of central galaxies by 
following the trajectories of individual galaxies in a fixed 
cluster potential. Approximately 30% of their simula- 
tions produced very massive galaxies at the cluster cen- 
ter, which they identified as cD galaxies. In their sim- 
ulations, which treated dynamical friction, mergers and 
tidal stripping simultaneously, the growth of the cD oc- 
curred throughout the lifetime of the cluster. Merritt 
(1984; 1985) reached a rather different conclusion using 
a statistical description of the galaxy orbital energy and 
mass distribution, and a Fokker-Planck treatment of the 
evolution of these distributions due to tidal stripping and 
dynamical friction. He found that neither tidal stripping 
nor mergers occur at a significant rate in rich clusters 
after virialization, and concluded that the morphology 



of BCGs in a rich cluster and of cDs in particular 
is fixed during the early stages of cluster evolution, 
probably during the collapse and virialization of com- 
pact groups or poor clusters. In the post-collapse stages 
of evolution, the central galaxy was found to grow in lu- 
minosity by only ~ IL*. Dynamical analysis of galaxies 
with multiple-nuclei led Merrifield & Kent (1991) and 
Blakeslee & Tonry (1992) to conclude such galaxies have 
probably increased their luminosities by no more than « 
1-2L* over a typical cluster lifetime as a result of galactic 
cannibalism. These findings are consistent with the mod- 
est post-virialization growth rates advocated by Merritt 
(1985). 

In principle, the formation of central galaxies in clus- 
ter and group environments is best studied with high- 
resolution, N-body simulations {e.g., Funato, Makino & 
Ebisuzaki 1993: Bode et al. 1994; Garijo, Athanassoula 
& Garcfa-Gomez 1997). The results of such simulations, 
however, depend sensitively on the assumed initial con- 
ditions, so it is important to start with a physically- 
motivated paradigm of structure formation. The first 
attempt to do so was by Dubinski (1998), who followed 
the formation and evolution of the central BCG within 
a single rich cluster chosen from a cold-dark-matter sim- 
ulation. In his simulation, which included both stellar 
and dark matter components, a massive central central 
galaxy formed via the merger of several massive "proto- 
galaxies" early in the cluster history. Due to the filamen- 
tary structure of the cluster at these times, the mergers 
themselves were found to be highly anisotropic, with the 
protogalaxies infalling radially along with filaments. In 
agreement with previous findings, Dubinski (1998) found 
that cannibalism by dynamical friction during the post- 
collapse evolution accounts for only a small fraction of 
the accreted mass. It is worth noting, however, that the 
central galaxy in his simulation did not develop the ex- 
tended envelope that is characteristic of cD galaxies. 

Another, quite different, formation model for cD galax- 
ies involves cooling flows (see, e.g., Fabian 1994). Since 
the radiative cooling times for intracluster gas are short 
enough that gas can cool and settle to the cluster center, 
it has been suggested that cD envelopes may arise from 
the gradual deposition of this cool gas. However, this 
possibility now seems remote in light of XMAI-Ncwton 
observations that the X-ray gas does not cool signifi- 
cantly below a threshold temperature of kT ^ 1-2 keV 
(Kaastra et al. 2001; Peterson et al. 2001; Tamura et al. 
2001). 

Clearly, an understanding of the mass distribution of 
the galaxies and protogalaxies from which cD galaxies 
were assembled would be helpful in discriminating be- 
tween the various models. Globular clusters (GCs) are 
useful in this regard since they are bright (with a mean 
absolute magnitude of My — —7.3 mag; Harris 2001) 
and are found in large numbers around massive, early- 
type galaxies. Moreover, observations of nearby dwarfs 
and giant ellipticals have shown that there is a tight 
correlation between mean GC metallicity and the total 
galaxy luminosity (e.g., Brodie & Huchra 1991). Thus, 
GCs should be powerful diagnostics of the cD formation 
process since the total luminosities and masses of the 
cannabilized galaxies (assuming cD galaxies formed in 
this way) must be imprinted in the metallicities of their 
present-day GC systems. Yet color (and thus metallic- 
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ity) information is available in the literaure for the GC 
systems of only five cD galaxies''. 

In this paper, we present B and / photometry of the 
GC systems of four distant cD galaxies and place con- 
straints on the distribution of progenitor masses. This 
sample nearly doubles the number of cD galaxies with 
GC metallicity distributions. The paper is organized as 
follows. In §2, we describe the observations and data 
reduction procedures. §3 presents a summary of the ob- 
served properties of the galaxies and their GC systems. 
In §4, we present an analysis of the observed GC metal- 
licity distribution functions. We summarize our findings 
and draw conclusions in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

We used WFPC2 on board HST to obtain F450W 
(broadband B) and F814W (broadband /) images of 
four cD galaxies: NGC 541, NGC 2832, NGC 4839 and 
NGC 7768. The observations, corresponding to program 
GO-8184, are summarized in Tabled Three of the galax- 
ies (NGC 2832, NGC 4839 and NGC 7768) were selected 
from the sample of Schombert (1986). The fourth galaxy, 
NGC 541, was selected from the compilation of Bird 
(1994), although it has in the past been variously classi- 
fied as D or cD {e.g., Dressier 1980; Bird 1994). To clarify 
the issue of its morphological type, we subsequently ob- 
tained B and / images of NGC 541 with the VLT. These 
data and our analysis of them are presented in i)2.1.1l 
Properties of the host clusters for each of our program 
galaxies are summarized in Table 

The raw WFPC2 data were processed with the stan- 
dard STScI pipeline using the best available calibration 
files. We verified that the images were properly aligned 
with cross-correlation measurements of the shifts using 
the DITHER package in IRAF^ Cosmic rays were re- 
jected by combining images of a given galaxy and filter 
using the IRAF task GCOMBINE. 

2.1. Surface Photometry 

The images were "mosaiced" using the task WMO- 
SAIC in the STSDAS package. Bad pixels, regions af- 
fected by vignetting, and any remaining cosmic rays were 
masked. The images were visually inspected and any 
bright contaminating objects present in the fields were 
also masked. Surface photometry was performed with 
the task ELLIPSE, which is based on an algorithm by 
Jedrzejewski (1987). The algorithm was run on the / 
image and the resulting solution was transferred to the 
B band image, where ELLIPSE was run in "no-fit" mode. 
The conversion to Johnson B and / magntiudes was done 
in the same manner as for the point sources (see ^2.2|l ex- 
cept that an additional 0.1 mag correction was applied 
to correct the zeropoints to infinite aperture (Holtzman 
et al. 1995). 

* Based on the compilation of van den Bergh (2001). Note that 
we include M87 as a cD galaxy (de Vaucouleurs &; Nieto 1978) and 
have added NGC 6166 (Bridges et al. 1996) to his list. 

^ Throughout this paper, we identify NGC 4938 with the sub- 
cluster projected ft: 1.1 Mpc to the southwest of the Coma cluster 
center. 

® IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



A common difficulty in performing surface photometry 
with WFPC2 is the limited field of view, which greatly 
complicates the estimation of the background sky bright- 
ness. Simply stated, the galaxies are still apparent even 
in the outskirts of our images. One approach to this 
problem is to fit a parametric curve (such as a de Vau- 
couleurs law or a power-law model) plus some constant 
background to a section of the profile, and estimate the 
sky directly from the fit. In our case, such a procedure 
is unsatisfactory. For the cD galaxies examined here, it 
is risky to assume a specific functional form when there 
is no a priori knowledge of where envelope will appear. 
Furthermore, sky measurements arrived at in this way 
will have independent errors in different bandpasses, po- 
tentially biasing the inferred color profiles at faint levels. 
Finally, as emphasized by Schombert (1986), no paramet- 
ric function can capture the peculiarities of the profiles 
that may affect the sky estimates. 

In an attempt to alleviate these problems, we adopted 
the following procedure. In the outskirts of our fields, the 
galaxy contribution to the total intensity, though non- 
negligible, is nevertheless small. By estimating the ratio 
of the sky intensities directly from the ratio image of the 
B and I images, it is possible to reduce the uncertainty in 
the estimated sky color. Let us denote the counts due to 
the sky in the / and B bandpasses by Sj and Sb , respec- 
tively, and the corresponding intensities from the galaxy 
light by Gi and Gb- Quite generally, the galaxy intensi- 
ties in the different bandpasses will be certain fractions of 
the sky intensities: i.e., Gj = aSi and Gb = PSb- For a 
typical giant eUiptical, {B — I)gai ~ 2.2 [e.g., Fukugita, 
Shimasaku & Ichikawa 1995; see below also) whereas the 
color of the "sky" seen by HST is {B - I)sky ~ 1-5 [HST 
Data Handbook). In general, we then have 



log(<7) EE log(/3/a) = 0.4 



[B — I)sky — {B — I)gal 



(2) 



and, if we estimate the ratio, R, of the sky intensities in 
the outskirts of the fields. 



R 



Si + aSi Si f 1 + a 



Sb + PSb SB\^ + qoi 



(3) 



The value of R was calculated for a range in plausible 
values of a and q. The distribution of R is sharply peaked 
at i? ~ 1.023S'//S'b, with an interquartile range of ~ 0.01 
and a median of i? = 1.0255//S's. We therefore take 
i?/1.025 for our estimate of Si/Sb, with R estimated 
directly from the ratio image. 

Since the galaxies are fainter with respect to the sky in 
the bluer bandpass, Sb was determined from the outer 
regions of the B images, along the direction where the 
profile declines most steeply. The sky brightness in / was 
then obtained from Si = SbR/^.025. Note that in the 
regions where S b was measured, the galaxy contribution 
is at least a factor of 20 less than that of the sky. We fur- 
ther restricted the final surface brightness profiles to be 
no more than 2.5 mag fainter than the estimated sky to 
guard against an artificial drop in the profile due to back- 
ground contamination by galaxy light. This procedure 
was used for all galaxies except NGC 2832, which ex- 
tends well beyond the WFPC2 field even in the B band. 
In this case, the estimate of Sb would be badly biased so 
we fixed the _B-band sky value by matching our photom- 
etry to the S-band profile of Peletier et al. (1990). The 
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main advantage of this sky estimation method is that it 
guarantees a priori that the measured sky color will be 
close to its actual value. 

Surface brightness measurements for the four galaxies 
are presented in Table 13 along with the results of our 
isophotal analysis. Surface brightness profiles are shown 
in the leftmost panels of Figures ESI Profiles from the 
literature are shown for comparison. We stress that, in 
the case of NGC 2832, the agreement at faint levels is 
forced upon our photometry for the purpose of sky de- 
termination. For the remaining galaxies, the profiles are 
fully independent from those in the literature. In each 
case, the agreement is good except for the outer parts of 
NGC 541, where the profile of de Juan, Colina, & Perez- 
Fournon (1994) turns downward compared to ours. We 
speculate that this downturn might be a consequence of 
their overestimating the background level: since our sky 
estimates are taken directly from the WFPC2 images, 
they are unlikely to be underestimated. 

Radial profiles of ellipticity, major-axis position angle 
and the parameter^ are shown in the right panels of 
Figures^^ As is common among BCGs, the ellipticities 
of NGC 541, NGC 4839 and NGC 7768 tend to increase 
with galactocentric radius {e.g., Porter et al. 1991) while 
their major-axis position angles tend to align with those 
of the surrounding clusters (e.g., Plionis 1994; West 1994; 
Kim et al. 2002). The ellipticity of NGC 2832 also in- 
creases with galactocentric radius, but in this case, the 
galaxy's major-axis at large radii is oriented almost per- 
pendicularly to the surrounding cluster. Note that the 
profiles of NGC 541 and NGC 7768 flatten abruptly at 
small radii (r <^ 1"). As Figure [S] clearly shows, the flat- 
tening in both case is due to the presence of compact 
dust disks in their nuclei. Tran et al. (2001) find such 
compact disks in ~ 18% of the objects in their WFPC2 
survey of early-type galaxies. It is curious that two of 
the four cD galaxies in our sample contain such disks. 

2.1.1. NGC 54I: Comparison with VLT ground based 
photometry 

As mentioned in the cD classification requires the 
presence of an envelope that departs from an r^^'^-law 
at low surface brightness. Our HST profile for NGC 541 
shows such an envelope starting at ^ 10", but we note 
that this galaxy has in the past been classified as both 
cD (Bird 1994) and D (Dressier 1980). To better char- 
acterize the surface brightness profile at large radii, we 
obtained B and / images of NGC 541 with FORSl on 
the VLT on 12 December, 2001 (see Figure El . Expo- 
sure times were 180 and 20 seconds in J, and 420 and 
60 seconds in B. The short exposures were used to re- 
place those pixels close to the galaxy center that were 
saturated on the long exposures. The FWHM of iso- 
lated stars on the images was measured to be ~ 1.1" 
in B and ~ 1.4" in /. The ground-based data, which 
have much greater areal coverage than the HST data 
and thus offer a more direct means of measuring the lo- 
cal background level, were reduced using standard IRAF 
routines and surface photometry was performed with the 
task ELLIPSE. Because the observations were taken in 
non-photometric conditions, the HST photometry was 

^ Bit indicates "boxy" isophotes if negative and "disky" ones if 
positive (Jedrzejewski 1987). 



used to set the photometric zeropoints by matching the 
profiles in the range 4" <; r ^ 10". 

The HST and VLT profiles are compared in Figure 
An envelope is clearly seen in the ground-based pro- 
file, confirming the classification of NGC 541 as a cD 
galaxy. There is generally good agreement between the 
two datasets, except at the faintest levels where the HST 
profile shows a steeper decline. The increase in elliptic- 
ity for r ^ r^ that was suggested by the HST data in 
Figure n] is readily apparent in the VLT data. It is also 
apparent that, in the transition region between the PC 
and WF chips (typically around R ~ 25"), the isophotal 
parameters derived from the HST data are somewhat less 
reliable because of the limited and varying areal coverage 
in this region. 

As the ground-based surface brightness profile becomes 
shallower at large radii, i?4 increases dramatically. This 
behavior can be understood by inspecting Figure Be- 
yond r ~ 1', there is significant contamination by close 
companions, both to the northeast and southwest. In 
fact, the various profiles overlap, producing a rather lin- 
ear structure that runs diagonally through the field. As 
a result, the isophotes of NGC 541 have elongated ends 
that are aligned with this feature. This behavior, which 
is manifested as an increase in the B4 parameter for the 
outermost isophotes, is unlikely to be is an artifact the 
sky estimation, since the slope of the surface brightness 
profile changes rather abruptly, and setting the sky to 
obvious overestimates does not make the increase in B4 
disappear. Given the presence of significant structure 
in the background light, what we are seeing is probably 
a consequence of our using a single value to character- 
ize a varying background in the vicinity of NGC 541. 
The surface brightness profiles shown in Figure [T] should 
therefore be treated with some caution beyond r ~ 1'. 

2.2. Point Sources: Detection, Classification and 
Photometry 

GCs at the distances of our program galaxies appear as 
unresolved point sources, even in the PC frames, so their 
analysis is readily amenable to PSF-fitting photometry. 
The images were first corrected for geometric distortion 
(Holtzman et al. 1995), and detection frames were cre- 
ated by applying a ring median filter (Seeker 1995) with a 
radius of five pixels for the PC and four pixels for the WF 
frames in order to eliminate the underlying galaxy light 
and to detect those sources with scale-lengths smaller 
than the filter radius. Detections were performed in- 
dependently in both filters using SExtractor (Berlin & 
Arnouts 1996) after setting the detection threshold to 
three connected pixels with counts 2a above the local 
background in the detection image (which was convolved 
with a gaussian filter). For some of our frames, there was 
an insufficient number of suitable stars to build an PSF 
directly, so we created empirical PSFs from archival ob- 
servations and used these PSFs in the analysis of our 
program images.^ Approximately one hundred stars per 
filter and chip were used to construct quadratically vary- 
ing PSFs with radii of ten pixels for the PC and five 
pixels for the WF chips. Using these PSFs, photometry 

* The F814W and F450W PSF were constructed from images 
of the Sagittarius dwarf galaxy (GO6701), and Andromeda VI 
(G08272), respectively 
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was performed with DAOPHOT II (Stetson 1987; 1993), 
and the object Hsts in both the F814W and F450W fil- 
ters were matched. Only those sources detected in both 
filters were retained. 

Corrections for charge transfer efficiency were carried 
out using the prescriptions of Dolphin (2000), and stan- 



dard HST magnitudes were measured using the zero- 
points given in Holtzman et al. (1995). The conversion 
to Johnson B and / was done using the following trans- 
formations (J. Holtzman, private communication), where 

(B-/)i.3 = (S-/)-1.3, 



/ = 

B = 



0.009(5 - /)2 3 



\niFsi4W - 0.012(5 - /)i.3 - 0.003(5 - /)? 3 - 

(mF450W + 0.087(5 - 1)1.3 - 0.008(5 -1)1 3 
1^^450^ + 0.193(5 - /)i.3 - 0.047(5 - /)? 3 

I 



0.021 if (5-/)i.3<0 

0.021 if (5-/)i.3>0 

-0.126 if (5-/)i.3<0 

-0.126 if (5-/)i.3>0 



A reddening correction was applied using the values of 
E{B — V) obtained from the DIRBE maps of Schlegel, 
Finkbeiner & Davis (1998) and the extinction curve of 
CardeUi, Clayton & Mathis (1989). Additionally, a k- 
correction was applied to the colors as described in ti2.3l 
Since at this stage the candidate lists still contain signif- 
icant contamination from compact background galaxies 
and foreground stars, additional selection criteria were 
imposed to obtain "cleaned" lists of GCs. Specifically, 
we discarded objects whose DAOPHOT x statistic^ ex- 
ceeded X = 2 in both 5 and /, and whose SHARP param- 
eters^o fell outside the range -0.35 ^ SHARP ^ 0.35. 
We also discarded objects with magnitudes more than 
3(7 brighter than the peak of the GC luminosity function 
(GCLF) in either 5 or I, where cr ^ 1.4 mag is the typ- 
ical dispersion of the GCLF (Harris 2001). In doing so, 
we adopted Mb — —6.73 and Mj = —8.39 mag for the 
GCLF "turnover" magnitudes, displaced to a distance 
obtained from the mean cluster redshift using Hq = 72 
km s^i Mpc^i (Freedman et al. 2001). 

Artificial star tests were performed to calculate the 
completeness function and to assess the extent to which 
our culling procedures would reject bonafide GCs. A 
total of 10,000 stars were added to each chip, 50 stars 
each time so as not to alter significantly the crowding 
characteristics of our frames. Artificial stars had 5 mag- 
nitudes in the range 24 ^ 5 ^ 27 mag and B — I was 
randomly assigned a value assuming the color distribu- 
tion is a Gaussian with the observed mean and dispersion 
for each galaxy. The artificial frames were then subjected 
to the same reduction procedure as the real data. Based 
on these tests, we conclude that the above selection cri- 
teria reject < 3% of the artificial stars. Photometry for 
GC candidates in our four galaxies that have photomet- 
ric uncertainties less than 0.3 mag in both 5 and / is 
presented in Tables 0- [3 



2.3. k-corrections 

® The DAOPHOT x is a robust estimate of the ratio of the 
observed to expected scatter about the model profile (Stetson & 
Harris 1988). 

The DAOPHOT SHARP parameter is an image radius index 
which is greater than zero if the object is more extended than 
the expected stellar profile and less than zero when the detection 
appears sharper. The expected value for a point source is zero 
(Stetson & Harris 1988). 



Our program galaxies have redshifts in the range 
0.018 ^ z 0.027. As we now show, redshift ef- 
fects on the observed magnitudes and colors become sig- 
nificant at these distances. Since these effects depend 
rather strongly on the object spectrum, we derive below 
k-corrections for both the GCs and the host galaxies. 
In doing so, we make the assumption that the GCs are 
old objects, meaning that they are roughly coeval with 
GCs in the Milky Way. This assumption is supported by 
most observations of GCs belonging to luminous early- 
type galaxies in cluster environments (see Jordan et al. 
2002 and references therein). Their broadband colors are 
then determined mainly by metallicity, and not by age 
(Worthey 1994). 

The k-correction is defined as 

E{X,to)Si\)dX 



(4) 

where, following the notation of Poggianti (1997), 
E{X,to) is the luminosity measured at wavelength A 
at the present time to in the source rest frame, and 
S{X) is the filter transmission curve. In calculating 
k-corrections for the GCs, we used four template GC 
spectra, based on observations of Galactic GCs, hav- 
ing different metallicities and different levels of ultravi- 
olet emission^^ (Bonatto, Bica & Alloin 1995). Specif- 
ically, we used their types G2b/r, G3b/r, G4b/r, G5, 
which have respective mean metallicities of [Fe/H] — 
—0.4, —1.0, —1.5, —1.9 dex. To calculate the correction 
needed for solar-metallicity GCs, and for the galax- 
ies themselves, we used the elliptical template of Kin- 
ney et al. (1996), which is based on spectroscopy of 
NGC 1399, NGC 1404, NGC 6868 and NGC 7196. The 
resulting k-corrections are presented in Table |H1 

Figure IHl shows k-corrections, K{B — I)o, plotted as 
a function of rest frame (5 — /)o. The upper and 
lower sequences show the corrections for NGC 541 and 
NGC 7768, the galaxies with the lowest and highest red- 
shifts in our sample. It is clear that the variation with 
(5 — /)o is significant and must be accounted for. We 
take the corrected colors to be linearly related to the 
observed (5 — /)o by 

(5-/)o,k = {B-I)o-K{B~I)o = Ck{B-I)o+Dk (5) 

The results for cluster types with different ultraviolet behavior 
are very similar so the corrections were averaged. 
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where the constants are given in Table |5| 

In principle, we might also wish to correct the ob- 
served colors for evolutionary effects. However, any cor- 
rection for the fading stellar populations in the color of 
these objects will be small compared to the k-correction. 
For instance, using the populations synthesis models of 
Maraston (1998), we find that the expected evolution in 
{B - I)o between 13 and 14 Gyr is 6{B - I)q < 0.03 
mag for all metallicities. As the lookback time to our 
targets is < 0.3 Gyr, the evolution in color is expected 
to be S{B — /) < 0.01 mag — significantly smaller than 
the k-correction on the color. Note that the evolutionary 
effects on the observed magnitudes are, in general, not 
negligible compared to the k-corrections for our targets. 
However, there is a partial cancellation in the corrections, 
as both B and / fade with age (albeit at different rates). 
For the present work, we have applied a k-correction to 
just the colors, since we use the magnitudes only when 
calculating the specific frequencies in ii3.3l and then only 
the / band where the k-correction and evolutionary cor- 
rection go in opposite directions. Thus, our results for 
the GC specific frequencies of these galaxies are quite in- 
sensitive to the correction (i.e., the /-band magnitudes 
vary by < 0.03 over the redshift range of our program 
objects). 

3. ANALYSIS 

3.1. Color Magnitude Diagrams and Color Distributions 

A remarkable finding from the past decade of GC 
research is that a majority of giant elliptical galaxies 
have GC systems with bimodal color distribution func- 
tions {e.g., Gebhardt & Kissler-Patig 1999; Larsen et al. 
2001). Because the colors of old stellar populations are 
more sensitive to metallicity than age, this bimodality 
implies the presence of two subpopulations of differing 
metallicity. In this section, we examine the color distri- 
butions of the GC systems associated with our cD galax- 
ies, thus probing the GC chemical enrichment histories of 
early-type galaxies in their most extreme manifestation 
of luminosity and size. 

Figure El presents color magnitude diagrams for the 
four GC systems. All objects classified as unresolved by 
our culling procedure are plotted in this figure. Magni- 
tudes and colors have been corrected for reddening and 
extinction, and k-corrections have been applied to the 
observed colors. It is apparent that none of the GC sys- 
tems show obvious sequences arising from the presence of 
two subpopulations. This impression is confirmed in Fig- 
ure^] where we plot color distributions for all objects 
with (j{B — I) < 0.3 mag. The solid line is a kernel den- 
sity estimate of the color distribution (Silverman 1986) 
obtained with a normal kernel. The adopted smoothing 
parameter h was chosen using the prescription of Sil- 
verman (1986) which is well suited to a wide range in 
densities and is trivial to evaluate: h = 0.9An~^^^ where 
A — mm(standard deviation, interquartile range/1.34). 
The dashed curves show 90% confidence bands obtained 
with a smoothed bootstrap procedure (Silverman 1986). 

We have computed a variety of statistical indices for 
each of the color distributions shown in Figure lTHI includ- 
ing biweight estimates of location, Cbi, and scale, Sbi, 
Finch's (1977) asymetry index, AI, (akin to skewness) 
and Hoaglin, Mosteller & Tukey (1983) tail index, TI, 



(akin to kurtosis).^^ The results are summarized in Ta- 
ble where the upper and lower limits on each index 
represent 68% confidence bands obtained with a boot- 
strap procedure. To within the uncertainties, the width 
of the distributions is the same for all galaxies, and it 
is clear that none of the distributions show statisitically 
significant skewness or extended tails beyond what would 
be expected from a parent Gaussian distribution. Simi- 
larly, the DIP test (see Gebhardt & Kissler-Patig 1999) 
accepts the hypothesis of unimodality for each of the dis- 
tributions shown in Fiaure lTUI We caution, however, that 
although bimodality is formally rejected in each case, the 
photometric errors, a{B — I) ^ 0.2 mag, are sufficiently 
large that intrinsic bimodality (with an expected sepa- 
ration between the modes of 0.3—0.4 mag in _B — /; e.g., 
Forbes et al. 1998) would be difficult to detect given the 
current samples. 

3.2. Color Profiles and Ojfsets 

Mean colors for our program galaxies and their GC 
systems are given in Table 1111 which summarizes vari- 
ous observed and derived photometric data. From left 
to right, this table gives the galaxy name, adopted red- 
dening, average colors for the galaxy and GCs, the total 
y-band magnitude from RC3, the galaxy redshift, its 
total t^-band absolute magnitude, and approximate lu- 
minosity in units of i*. In Figure [TTl we plot the colors 
of the GC candidates as a function of galactoccntric dis- 
tance for each of the four galaxies. The solid line in 
each panel is the best-fit linear relation for GCs with 
a{B — I) < 0.25 mag. The dashed-lines are profiles of 
the galaxy light measured from our HST images. In all 
four cases, as is typical for GC systems, there is a nega- 
tive color gradient in the sense that the GC system be- 
comes bluer with increasing galactocentric distance {e.g., 
Strom etal. 1981). The GC color gradients parallel those 
of the underlying galaxies, with mean offsets ranging 
from 0.33 mag for NGC 541 to 0.50 mag for NGC 7768. 
If we make the somewhat daring assumption that the 
composite stellar populations of the galaxies follow the 
same color-metallicity relation as the globulars (Barmby 
et al. 2000), then these differences correspond to metal- 
hcity offsets in the range 0.8 A[Fe/H] 1.3 dex. 

Indeed, the color /metallicity offset apparent in Fig- 
ureEJis a ubiquitous feature that is observed in galaxies 
ranging from faint, low-surface dSph/dE galaxies to the 
brightest and most massive early-type galaxies. The left 
panel of Figure El shows mean GC and stellar metallic- 
ity plotted as a function of host galaxy magnitude for a 
large sample of early-type galaxies. Note that this plot is 
an updated version of Figure 1 of Cote et al. (2000), to 
which we have added elliptical galaxies from the sample 
from Kundu & Whitmore (2001) and our four cD galaxies 
(shown as the circled points). Improved stellar metallic- 
ity determinations for several dwarf galaxies have been 
taken from Shetrone, Cote & Sargent (2001) and Tol- 
stoy et al. (2003). For those galaxies from the survey of 
Kundu & Whitmore (2001), GC metallicities were calcu- 
lated by combining their {V — I) colors with the {V — /)- 
metallicity relation of Barmby et al. (2000). Metallicities 

For a description of the asymetry and tail indices and some 
of their properties, see Bird and Beers (1993). See Beers, Flynn & 
Gebhardt (1990) for a description of the biweight estimators. 
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for the galaxies themselves come from the integrated- 
light spectroscopy of Trager et al. (2000) when avail- 
able; for those galaxies lacking spectroscopic data, we 
estimated crude metallicities by combining their (B — V) 
colors from RC3 (de Vaucoulcm-s et al. 1995) with the 
{B — F)-metallicity relation of Barmby et al. (2000). 

It is clear from the left panel of Figure El that the 
GCs and stars obey reasonably tight relations that run 
roughly parallel to each other. The right panel of this 
figure shows, for those galaxies that have metallicity de- 
terminations for both the stars and the GC systems, the 
offset between the two components. Remarkably, this 
offset is nearly constant over the whole range of galaxy 
magnitude. An i^-test indicates that the offset appears 
to be constant over the full range in galaxy magnitude, 
with a weighted mean value of A[Fe/H]= 0.79±0.04 dex. 
We conclude that the metallicity offset between stars and 
GCs in these luminous cD galaxies is indistinguishable 
from that observed in normal ellipticals. This metallicity 
offset presents another apparently "universal" property 
of GC systems — analogous to the near-Gaussian nature 
of the GC luminosity function (Harris 2001) and the ap- 
parently constant efficiency of GC formation (Blakeslee 
et al. 1997; McLaughhn 1999a) — that must be ex- 
plained by any viable model of GC formation. 

3.3. Specific Frequencies 

The GC specific frequency, Sn = A^gc10°-^(*''^^+^^', is 
a convenient measure of GC formation efficiency relative 
to field stars in the underlying galaxy (Harris & van den 
Bergh 1981; cf. McLaughlin 1999a). Possible dependen- 
cies of Sn on morphological types have played an impor- 
tant role in the debate over how, and when, early-type 
galaxies form (Schweizer 1987; Ashman & Zepf 1992; van 
den Bergh 1995). Since the first "high-^Ar" giant galax- 
ies to be identified (M87, NGC 1399, NGC 3311) were aU 
cDs, it was believed that the cD environment was some- 
how conducive to the efficient formation of GCs relative 
to field stars (Harris 1988). However, with the acquisi- 
tion of more and better data, it became apparent that 
not all cD galaxies had high specific frequencies {e.g., 
Harris, Pritchet & McClure 1995), leading West et al. 
(1995) to propose a phenomenological model in which 
the GC specific frequency of BCGs correlates with the 
mass of the surrounding galaxy cluster. This suggestion 
was later confirmed observationally as a correlation be- 
tween specific frequency and cluster velocity dispersion 
or X-ray temperature (Blakeslee, Tonry & Metzger 1997; 
Blakeslee 1999). In other words, it appears that the prop- 
erties of the surrounding galaxy cluster are connected to 
both the cD envelope and the overall number of GCs, 
and that the mere presence of an envelope is not enough 
to guarantee a high GC specific frequency. 

To measure GC specific frequencies for our program 
galaxies, we must first estimate the sizes of their GC 
systems using the bright end of the GCLF, which we 
observe directly. For data such as ours, which do not 
reach the turnover, any attempt to fit simultaneously for 
the number of GCs, the turnover magnitude, and the 
dispersion of the GCLF would be confounded by large 
and highly correlated errors in the derived parameters 
(Hanes & Whittaker 1987; Seeker & Harris 1993). We 
therefore proceed by fixing the dispersion of the GCLF 
at a — 1.4 mag and adopting Mxo = —7.33 mag for the 



F-band absolute magnitude its turnover (Harris 2001). 

To obtain the completeness-corrected luminosity func- 
tion, we divide our GC sample into three radial regions: 
r < 20", 20" < r < 50" and r > 50" . We determine 
the completeness function, /(m), in each radial bin by 
fitting the function 

„, , 1 a(m — miirn) , n 

/(to) = - 1 V ^^ra^ (g) 

2 [ \/l + a^{m - mumY 

to the results of our artificial star tests, where mura is 
the magnitude at which / = 0.5 (Fleming et al. 1995). 
These functions were then used to correct the luminosity 
functions in each bin up to a magnitude obtained by 
averaging mura for the three regions. The three LFs 
were then added to find the completeness corrected LF 
for the full field. 
We fit for TVgc in a GCLF of the form 

N(m) = ArGc(27ra2)-i/2 exp[-(r7i - m^of I {la'^)] (7) 

where toto — {m — M) + A/to- Here Mto is absolute 
magnitude of the GCLF turnover and {m — M) is the dis- 
tance modulus obtained from the observed redshift (see 
Table El and i/o = 72 ± 8 km s^^ Mpc^^ (Freedman 
et al. 2001), which is consistent with the value obtained 
from WMAP (Spergel et al. 2003). We assume the un- 
certainty in the turnover magnitude is due to the un- 
certainty of 0.05 quoted in Harris (2001) plus the joint 
uncertainties arising from the distance determination (as 
a result of deviations from a pure Hubble fiow) and the 
uncertainty in the Hubble constant itself. The combined 
error in the turnover magnitude due to the latter effects 
is 

= (5//n(10))2(a,Jcz)2 + {h / H^ln{lQ)f (8) 

where we have adopted a constant value of Ucz — 300 km 
s~^ for the rms cluster peculiar velocity (Bahcall & Oh 
1996). Clearly, these uncertainties should be viewed as 
lower limits since background contamination, photomet- 
ric errors for the GCs, and model assumptions have not 
been included. 

We take A^qc in equation (7) as our best estimate for 
the total number of GCs within our field. Since this is 
a metric quantity, it is denoted as Nqq hereafter. For 
each galaxy, S'J^°* is measured by calculating the inte- 
grated magnitude of the underlying galaxy within our 
WFPC2 field with the aid of the B- and /-band galaxy 
models presented in § 12.11 To convert these B- and /- 
band magnitudes into the V bandpass, we assume typi- 
cal colors for ellipticals as given in Fukugita et al (1995): 
{B-V) = 0.96 mag and {V - I) = 1.31 mag. Table [HI 
summarizes the GC specific frequency measurements for 
our cD galaxies. From left to right, this table records the 
galaxy name, absolute magnitude, total number of glob- 
ular clusters, and GC specific frequency. Once again, 
these are all metric quantities. 

It is interesting to note that none of these galaxies 
show high specific frequencies, their values being typi- 
cal for ellipticals and even somewhat low for NGC 7768. 
Blakeslee et al. (1997) measured specific frequencies for 
NGC 2832, NGC 4839 and NGC 7768 as part of their 
survey of BCGs. Within the errors, our results are con- 
sistent with theirs for all three objects. NGC 7768 was 
also studied by Harris, Pritchet & McClure (1995), who 
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obtained Sn = 3 ± 2. Their result scales with their 
assumed value of Sn = 15 for M87. In light of more 
recent measurements the GCLF parameters in M87 by 
Whitmore et al. (1995), we adopt Sn = 11 for this 
galaxy (using the correction factor presented in equation 
17 of McLaughlin, Harris & Hanes 1994). The Harris 
et al. (1995) measurement for NGC 7768 then falls to 
Sn ~ 2.3. This is in good agreement with our determi- 
nation of Sn = 2.2±0.8. Finally, Marm-Franch & Apari- 
cio (2002) find Sn = 7.0 ± 1.9 for NGC 4839 — shghtly 
higher than, but still consistent with, our measurement 
of Sn = 5.6 ± 1.9. This is the first determination of S'jv 
for NGC 541 and its value — like those of the two other 
galaxies for which X-ray luminosities of the host clus- 
ter are available (see Table (SJ — is consistent with that 
expected from the cluster's X-ray temperature and the 
Sn-Lx correlation presented by Blakeslee et al. (1997). 

4. CONSTRAINTS ON THE MERGER HISTORIES OF CD 
GALAXIES 

As discussed in the physical properties of cD galax- 
ies, coupled with their unique locations near the dynam- 
ical centers of rich galaxy clusters, suggest that some as- 
pects of their formation and/or evolutionary histories dif- 
fer from those of normal elliptical galaxies. A wide range 
of evolutionary phenomena has been invoked to explain 
their properties, including cannibalism, tidal stripping, 
and star formation in coolings flows. Clearly, a resolution 
of the debate concerning which, if any, of these processes 
has played a role in the formation of cD galaxies would 
benefit from a knowledge of their merger histories. 

4.1. Protogalactic Mass/Luminosity Spectra From 
Globular Cluster Metallicity Distributions 

A method of using GC metallicity distributions to de- 
termine protogalactic mass spectra has been described by 
Cote et al. (1998; 2000; 2002). A fuU description of the 
method is presented in these papers, so here we simply 
note its basic assumptions and limitations. In brief, the 
method assumes that the protogalactic luminosity (and 
mass) spectrum can be approximated by a Schechter 
function (Schechter 1976), and that the present-day GC 
specific frequency of any surviving protogalactic "frag- 
ments" is Sn ~ 4. It is further assumed the galaxy under 
consideration has been assembled hierarchically, in a pri- 
marily dissipationless manner, so that no GCs are formed 
during mergers and interactions. In other words, star and 
GC formation during mergers are ignored entirely in this 
approach. The metallicities of individual GCs are instead 
taken to reflect differences in the depths of the gravi- 
tational potential wells in which they formed. Indeed, 
the tight empirical relation between mean GC metallic- 
ity and total galaxy magnitude, such as that shown in 
Figure 1121 is prima facie evidence for the importance of 
local environment in the chemical enrichment of GC sys- 
tems. As a consequence, we adopt a GC metallicity-host 
galaxy luminosity relation that is based on nearby dwarf 
galaxies and the metal-rich (bulge) GC populations of 
M31 and the Milky Way (Cote et al. 2002). 

The implementation of the model proceeds as follows. 
Using the above model inputs, simulated GC metallicity 
distributions are generated over a grid in two parameters 
chosen to describe the mass spectrum^'^ of protogalactic 

In converting to masses, we assume a constant mass-to-light 



fragments from which the galaxy was assembled: (1) the 
slope, a, of the Schechter function from which the sim- 
ulated protogalactic fragments were drawn; and (2) the 
ratio of the mass of largest protogalactic fragment, A4^, 
to the final mass, A^/ of the assembled galaxy: 

C^M^/Mf. (9) 

For each simulated GC system, we extract at random an 
identical number of GCs as contained in the actual sam- 
ple and add the appropriate amount of measurement er- 
ror to the simulated colors. The transformation between 
color and metallicity is performed using the transforma- 
tions of Barmby et al. (2000). We restrict the sample of 
GCs to lie on the range 1.23 < {B - /)o,k < 2.49 mag, 
which is equivalent to —2.5 < [Fe/H] < 1.0 dex. This 
generous range is expected to include the vast majority 
of bona fide GCs in these galaxies (see, e.g., Beasley et al. 
2000; Cohen, Blakeslee & Cote 2003). 

For each simulated GC system, a Kolmogorov-Smirnov 
(KS) test is performed to test the hypothesis that the 
simulated and observed metallicity distributions have 
been drawn from the same parent distribution. Simula- 
tions are carried out 100 times at each point in the a-C 
plane, giving a grand total of 60,000 simulations over the 
full grid. At each grid element, the average KS probabil- 
ity is recorded. The resulting probability surface, shown 
in Figures HSl CSl El and El gives the likelihood that 
the observed GC color/metallicity distributions can be 
reproduced with a protogalactic spectrum characterized 
by that a-^ pair. The four crosses in each figure mark 
representative simulations for each galaxy that have a 
better than 99% probability of being drawn from the 
same distribution as the observed GC metallicity distri- 
bution function. The upper panels of Figures 1141 1161 IT^ 
and 1201 show comparisons of the observed GC metallicity 
distribution functions with the results from these four 
simulations; the lower panels of these figures show the 
corresponding protogalactic mass spectra. 

The key result from these simulations is that the ob- 
served GC color/metallicity distributions seem to be per- 
fectly consistent with those expected if cD galaxies form 
through the accretion of numerous galaxies and proto- 
galactic fragments that managed to form the bulk of the 
GCs before being integrated into the central galaxy. At 
the same time, Figures [T!^20l reveal that a wide variety 
of protogalactic mass spectra are able to reproduce the 
observed metallicity distributions. For comparison, Cote 
et al. (2002) carried out an analysis of the 28 early-type 
galaxies of Kundu & Whitmore (2001) using the method 
described above, finding {() = 0.25 ± 0.03 (mean error) 
and a rather steep slope of (a) = —1.88 ± 0.03 (mean 
error). While these values are certainly consistent in all 
four cases with those found here. Figures I13H20I show 
that, with the possible exception of NGC 4839, shallower 
slopes are also allowed. Aside from NGC 4839 (which we 
discuss in detail below), the locus of most probable a-C 
pairs for each galaxy tend to form an elongated ridge that 
stretches from the lower-left corner of the plot {i.e., low ( 
and shallow slopes) to the middle-right region {i.e., inter- 
mediate C and steep slopes). The weaker constraints on 
a and C compared to those found by Cote et al. (2002), 
are a direct consequence of the fact that the GC metal- 
ratio of Tv = 5 for all protogalactic fragments. 
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licity distributions for these cD galaxies are dominated 
by photometric errors. With the current data, we are 
able to conclude that the GC metallicity distributions 
are consistent with the cDs being assembled via mergers, 
accretion and tidal stripping, but we are unable to place 
firm constraints on the detailed shape of the mass spectra 
of the cannibalized objects (with the possible exception 
of NGC 4839; see below). We now turn out attention to 
the mechanism(s) responsible for these putative mergers. 

4.2. Dynamical Friction and the Luminosity Function 
of Cannibalized Galaxies 

Note that above discussion refers implicitly to the ini- 
tial luminosity function of galaxies, protogalaxies and 
protogalactic "fragments". Since our simulations treat 
mergers in a purely statistical manner, equal merger 
probabilities are assumed for all galaxies and protogalac- 
tic fragments. This assumption is clearly an oversimpli- 
fication: if mergers proceed through dynamical friction, 
the merger probability will be a function of mass {e.g. 
Binney & Tremaine 1987). In this section, we examine 
how dynamical friction is expected to modify the initial 
luminosity function. Our aim is compare the luminosity 
function of cannibalized galaxies predicted by dynamical 
friction arguments with that obtained in our analysis of 
the GC metallicity distributions. 

At the centers of clusters where cD galaxies reside, the 
galaxy luminosity function will evolve with time because 
the rate of dynamical friction is mass (and, hence, lumi- 
nosity) dependent. If cannibalism as a result of dynami- 
cal friction is the primary mechanism by which cD galax- 
ies evolve, then we would expect the luminosity function 
of cannibalized galaxies to be skewed toward luminous 
galaxies. Note that, for our purposes, we are not con- 
cerned with the rate at which galaxies sink to the cluster 
center, but merely with the final luminosity distribution 
of all galaxies cannibalized over the cluster lifetime. 

We assume that the initial luminosity function is given 
by a Schecter function. 



(j)i{L) = n^{L/L^)°'^ exp(— L/L*)L^ 



(10) 



and that the joint luminosity and spatial density distri- 
bution of galaxies within the cluster, n{r, L) is given by 
the product of (j)i{L) with a spatial distribution ry(r), such 
that ri(r, L) — (j)i{L)'q[r) . A galaxy orbiting on a circular 
orbit in an isothermal potential reaches the cluster center 
after a time (Binney & Tremaine 1987) 



1.65rfa 
GA/lnA' 



(11) 



Here rt is the initial distance of the galaxy, a is the cluster 
velocity dispersion, and In A is the standard Coulomb 
logarithm. Setting tf to the age of the cluster, Tdus, we 
can solve for the distance, ri{M), within which galaxies 
of mass M will have spiraled to the center. Assuming a 
constant mass-to-light ratio, T, we can also express this 
distance as a function of L, 



Tr,,„,GlnA L 



n 0.5 



1.65 



(X L 



0.5 



(12) 



For the host clusters of our cD galaxies, r.i varies be- 
tween ~ 13 and ^ 1300 kpc for galaxy masses in the 



range 10^ to 10^^ AIq. The luminosity function of can- 
nibalized galaxies will be given by 



0,(i) =4^4(L) 



n(L) 



ri{r)r'^ dr. 



(13) 



The surface density profile of galaxies in the inner re- 
gions of clusters goes as p{r) oc r~^ (e.g., Beers & Tonry 
1986; Merritt & Tremblay 1994; McLaughhn 1999b). 
For a spherical cluster, this corresponds to a three- 
dimensional density distribution of r]{r) oc r^^, which 
in turn implies 



4>ciL) oc (t>i{L)L 



0.5 



(14) 



In other words, the slope of the luminosity function of 
cannibalized galaxies, ac, is related that of the initial 
luminosity function, a^, through the relation 



ai + 0.5. 



(15) 



This expression provides a a crude connection between 
the expected luminosity function of cannibalized galaxies 
to that of the surrounding cluster at the time of its for- 
mation under the simplifying assumptions made. Note 
that except in the case of a very steep initial density 
distribution, ac ^ ai regardless of the exact value of 
ac- ai. 

4.3. Comparing the Luminosity Functions 

Is it plausible that these cD galaxies were assembled 
mainly through dynamical friction? We may examine 
this possibility by assuming, for the time being, that the 
protogalactic luminosity functions found in M4.1I are the 
outcome of cannibalism through dynamical friction, so 
that a ~ Qic and ai = a — 0.5. 

For NGC 54 1, NGC 2832 and NGC 7768, Figures d 
[THI and I19I20I show that a rather wide range of proto- 
galactic luminosity functions are capable of reproduc- 
ing the observed GC metallicity distributions in these 
galaxies. Depending on the precises values of C, slopes 
in the range —2 < a < —0.75 are capable of yield- 
ing simulated metallicity distributions that are in close 
agreement with those observed; the constraints on a 
are somewhat tighter in the case of NGC 4839, with 
—2 <^ a <^ —1.1. These values of a would require the 
host clusters (i.e., Abell 194, Abell 779 and AbeU 2666) 
to have initial luminosity functions with slopes in the 
range —2.5 <^ ai <^ —1.25. For Abell 1656, the corre- 
sponding limits are —2.5 <^ ai <^ —1.65. Are such slopes 
feasible, or can they be ruled out from observations of 
the present-day luminosity functions in these clusters? 

There is, to the best of our knowledge, no published 
luminosity function for Abell 779, so we restrict ourselves 
to Abell 194, AbeU 1656 and Abell 2666. Interestingly, 
there have been recent reports of steep luminosity func- 
tions in both Abell 194 and Abell 2666. Trentham (1997) 
finds Qfobs > -1.6 in Abell 194, while de Propris et al. 

(1995) claim aobs 2.2 for Abell 2666. We caution, 

however, that in both cases, the measured luminosity 
functions are uncertain to the point being essentially un- 
constrained {e.g., Trentham quotes a best-fit value of 
Q^obs — —2.2^1^ for Abell 194). Firm conclusions on 
the viability of dynamical friction as the main mecha- 
nism responsible for the formation of these galaxies must 
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await the measurement of improved GC metallicity dis- 
tributions and, equally important, accurate luminosity 
functions for their host clusters. 

Observational constraints on the luminosity function in 
Abell 1656 are of much higher quality. Several compre- 
hensive studies of the luminosity function in this clus- 
ter which have yielded results that are in good agree- 
ment: aobs = -1.42 ± 0.05 (Bernstein et al. 1995), 
aobs = -1.41 ± 0.05 (Seeker, Harris & Plummer 1997), 
ttobs = —1.33 ±0.06 (Beijersbergen, Schaap & van der 
Hulst 2002), and ctobs = -1.31 ± 0.05 (Mobasher et al. 
2003). Strictly speaking, these measurements refer to 
the global luminosity function of the Abell 1656 clus- 
ter and not the subcluster that hosts NGC 4839, but 
Mobasher et al. (2003) find the luminosity function in 
the NGC 4839 subcluster to be indistinguishable from 
that of the cluster as a whole. The average of the above 
values, (aobs) = —1.37, falls outside the range of proba- 
ble tti values given above {i.e., —2.5 ^ at <^ —1.65). We 
tentatively conclude that the hypothesis that cD galax- 
ies were assembled mainly through mergers driven by 
dynamical friction seems inconsistent with the available 
data for NGC 4839 and Abell 1656, the galaxy/cluster 
pair in our sample with the highest quality GC metal- 
licity distribution and luminosity function. Recall from 
that various dynamical evidence suggests that cD 
galaxies typically grow by only « 1-2 L* as a result of 
cannibalism through dynamical friction during the post- 
virialization regime {e.g., Merritt 1985; Lauer 1985; Du- 
binski 1998). Since NGC 4839 has a luminosity of ^ 7L*, 
this conclusion is in agreement with the expectation from 
the dynamical evidence that dynamical friction alone can 
not account for the observed luminosity. 

5. CONCLUSIONS 

We have nearly doubled the number of cD galaxies with 
available color/metallicity distributions for their GC sys- 
tems, thus probing the realm of early- type galaxies in 
their most extreme manifestation of luminosity and size. 
Despite the undeniably special nature of these galaxies, 
their GC systems seem remarkably normal in all the 
properties we have explored: specific frequency, mean 
metallicity, and the metallicity offset between the GCs 
and host galaxy stars. An analysis of the protogalactic 
mass spectra for these galaxies, while more model depen- 
dent, similarly reveals no obvious anomalies compared to 
normal ellipticals. 

The GC specific frequencies of these cD galaxies fall 
within the expected range given the known scaling rela- 
tions between specific frequency and cluster X-ray tem- 
perature and velocity dispersion. Our specific frequen- 
cies corroborate previous measurements for three of our 
galaxies; the measured specific frequency for the remain- 
ing object, NGC 541, is consistent with expectations 
given the mass of Abell 194, the surrounding cluster. 
Thus, our findings are consistent with a picture in which 
the efficiency of GC formation scales with the total mass 
of the cluster itself, implying that the GC formation ef- 
ficiency per unit total mass {i.e., including dark matter) 
is universal (Blakeslee et al. 1997; Blakeslee 1999). It 
would be interesting to measure density profiles for the 
X-ray emitting gas these clusters, to see if the total num- 
bers of GCs in these galaxies are also consistent with 
a constant efficiency of GC formation per baryon mass, 



as suggested by McLaughhn (1999a). In any event, it 
seems clear that the GC systems belonging to these cD 
galaxies trace the cluster properties so that, as Blakeslee 
(1999) notes, it may be more correct to view "high-S'Ar" 
galaxies as being underluminous given their total masses, 
rather than ascribing to them an unusually high GC for- 
mation efficiency. Possible explanations for this "miss- 
ing light" problem include tidal heating of the gas as the 
galaxy cluster collapses (Blakeslee et al. 1997), and star- 
bust driven galactic winds (Harris, Harris & McLaughlin 
1998) which could be responsible for halting star for- 
mation in BCGs to an extent that is controlled by the 
cluster's mass. 

Our results suggest that, whatever mechanism is re- 
sponsible for the formation of cD galaxies, perhaps mak- 
ing some of them appear underluminous compared to 
their host clusters in the process, it seems to leave the 
GCs unscathed. Indeed, the metallicity offset between 
the CG system and the underlying galaxy is indistin- 
guishable from that observed for other early-type galax- 
ies, and thus presents us with another uniformity which 
must be tied in a fundamental way to their formation 
process. A plausible origin for this offset may rely on 
the timing of GC formation. We find the mean metallic- 
ity of the gas out of which GCs formed to be roughly six 
times lower than that of galaxy itself; if the GCs form be- 
fore the bulk of the stars, then a metallicity offset seems 
a natural, if not inevitable, outcome. This admittedly 
qualitative scenario agrees with the idea that star forma- 
tion was inhibited for some central galaxies but not their 
GC systems, as the GCs might have formed prior to the 
onset of this as-yet- unidentified process. 

Different models for the formation of cD galaxies dif- 
fer mainly on the timescales over which the cD galaxies 
form. We have examined the possibility that cD galaxies 
have been assembled through the accretion of numerous 
galaxies and protogalactic fragments that managed to 
form their GCs prior to being captured and disrupted 
by the central galaxy. Our key finding in this regard is 
that such a process appears to be entirely consistent with 
the observed GC metallcity distributions for these galax- 
ies. Unfortunately, a wide variety of protogalactic mass 
spectra are able to reproduce the observed metallicity 
distributions given the rather rather poorly constrained 
metallicity distributions. 

What mechanism might be responsible for such merg- 
ers? In the case of NGC 4839, the cD in the southwest 
extension of Abell 1656 and the galaxy in our sample 
having the most reliable GC metallicity distribution, we 
find the present-day luminosity function of Abell 1656 
to be too shallow to be consistent with the assembly 
of the galaxy from dynamical friction. In other words, 
while it is possible to choose an ensemble of galaxies, 
protogalaxies and protogalactic fragments that, once as- 
sembled, would produce a GC system with the appro- 
priate metallicity distribution, the required luminosity 
function of these objects seems overly steep compared to 
that measured at the present time. This suggests that, 
in agreement with the previous conclusions based on dy- 
namical arguments (Merritt 1985; Lauer 1985), cannibal- 
ism through dynamical friction is not the primary means 
by which cD galaxies are assembled. The properties of 
the GC system instead favor a scenario {e.g., Dubinski 
1990) in which the cDs form rapidly, and at early-times. 
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via hierarchical merging prior to cluster virialization. 
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Fig. 1. — (Left Panel) Surface brightness profiles for NGC 541 in B (filled squares) and / (filled circles). The squares are a Gunn-Thuan 
r-band profile from de Juan et al. (1994). The dot-dashed lines are best-fit de Vaucouleurs laws in the range 1.4 < ri/" < 1.8. Note the 
flattening inside r ~ 1" caused by the central dust disk (see Figure |SJ. (Right Panels) EUipticity (top), position angle (middle) and S4 
(bottom) as a function of galactocentric distance. The arrow in the middle panel indicates the position angle of the major axis of Abell 194 
(Lambas, Groth & Peebles 1988; Plionis 1994). 
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Fig. 2. — ('ie/f Panel) Surface brightness profiles for NGC 2832 in B (filled squares) and / (filled circles). The squares and circles 
are respective B- and ij-baiid profiles from Peletier et al. (1990). Crosses denote the photographic profile of Schombert (1986). The 
dot-dashed lines are best-fit de Vaucouleurs laws in the range 1.4 < r^/* < i g (Right Panels) Ellipticity {top), position angle {middle) 
and B4 {bottom) as a function of galactocentric distance. The arrow in the middle panel indicates the position angle of the major axis of 
Abell 779 (PUonis 1994). 
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Fk;. 3. (Left Panel) Surface brightness profiles for NGC 4839 in B (crosses) and / (filled circles). The squares and circles are respective 
B- and r-band profiles from J0rgensen, Franx & Kjaergaard (1992). Crosses denote the photographic profile of Schombert (1986). The 
dot-dashed lines are best-fit de Vaucoulcurs laws in the range 1.4 < r^/'^ < 1.8. (Right Panels) Ellipticity (top), position angle (middle) 
and B4 (bottom) as a function of galactocentric distance. The arrow in the middle panel indicates the position angle of the major axis of 
Abell 1656 (Rhee, van Haarlen & Katgert 1991; PUonis 1994). 
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Fig. 4. — (Left Panel) Surface brightness profiles for NGC 7768 in B (crosses) and I (filled circles). The squares are a V-band profile from 
Malumuth & Kirshner (1985), while the crosses denote the photographic profile of Schombert (1986). The dot-dashed lines are best-fit de 
Vaucouleurs laws in the range 1.4 < r^/^ i g Note the flattening inside r ~ 0''5 caused by the central dust disk (see Figure|3. (Right 
Panel) EUipticity (top), position angle (middle) and B4 (bottom) as a function of galactocentric distance. The arrow in the middle panel 
indicates the position angle of the major axis of Abell 2666 (Lambas, Groth & Peebles 1988). 
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Fig. 5. — F450W images of the central regions of NGC 541 and NGC 7768. Note the obvious dust disks. 




Fig. 6. — /-band image of NGC 541 taken with FORSl at the VLT. North is up and East is to the left in this image, which measures 
6.8' X 6.8'. Overlayed are the fitted eUipses for NGC 541, along with the WFPC2 field of view. 
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Fig. 7. — (Left Panels) Comparison of HST and VLT surface photometry for NGC 541 in the B (bottom) and / (top) bandpasses. The 
dotted line shows the sky brightness in the VLT frames; the dashed line is the corresponding value in the HST frames. The VLT profile 
was shifted to agree with that from HST in the region 1.4 < ri/4 < 1.8 (i.e., between 4" and 10"). The dot-dashed Une is the best-fit 

de Vaucoulcurs law in this range. The cD envelope is clearly visible as the excess above this line, beginning at fa 10". (Right Panels) 
Comparison of the ellipticity (top), position angle (middle) and B4 (bottom) measurements as a function of r'^/*. 
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Fig. 9. — Color magnitude diagrams for globular cluster candidates in our four program galaxies. The magnitudes and colors have 
been corrected for reddening and extinction; k-corrections have been applied to only the latter. The dashed lines show the color range, 
1.1 < (S — /)o < 2.4, used to isolate globular clusters from foreground stars and background galaxies. 
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program galaxies. The solid line is the line of best fit for objects with a{B — I) < 0.25 mag. The dashed line is the color profile of the 
underlying galaxy light. 
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Fig. 12. — (Left Panel) Mean metallicity of globular clusters (small filled circles) and stars (small open circles), plotted as a function 
of host galaxy absolute magnitude. Circled points indicate the four cD galaxies examined here. The dashed lines arc the lines of best-fit 
for the two samples. (Right Panel) Metallicity offset A[Fe/H] between stars and globular clusters, plotted as a function of host galaxy 
absolute magnitude. The dashed line indicates the weighted average: A[Fe/H] = 0.79±0.04 dex. Circled points indicate the four cD 
galaxies examined here. 
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J J 



IQB IQ9 IQIO 10" 1012 

Fig. 14. — (Upper Panel) Cumulative metallicity distribution of 208 globular cluster candida tes a ssociated with NGC541 (black curve). 
The four colored curves corresponds to the simulations marked by the colored crosses in Figure [T^ In each case, a KS test indicates that 
there is a better than 99% probability that it was drawn from the same parent distribution as the real data. (Lower Panel) Protogalactic 
mass spectra corresponding to the simulated metallicity distributions shown above. The dotted lines show powerlaw mass functions, 
A'^(A^) oc , with exponents of 7 = —2, —1.6, and —1.2. The vertical arrow denotes the total mass of the final galaxy. A constant 

mass-to-light ratio of Tv = 5 has been assumed in the conversion of luminosities to masses for both the final galaxy and the protogalactic 
fragments. 




Fig. 15.— 



Same as Figure El except for NGC 2832. 
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Fig. 16. — Same as Figure lT^ except for NGC 2832. The cumulative metallicity distribution is based on a total of 323 candidate globular 
clusters. 
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Fig. 17.— 



Same as Figure El except for NGC 4839. 
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Fig. 18. — Same as Figure lT^ except for NGC 4839. The cumulative metallicity distribution is based on total of 350 candidate globular 
clusters. 




Fig. 19.— Same as Figure [T^ except for NGC 7768. 
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Fig. 20. — Same as Figure lT^ except for NGC 7768. The cumulative metallicity distribution is based on total of 123 candidate globular 
clusters. 
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Table 1. Observing Log for GO 8184 



Galaxy 


R.A. 


Dccl. 


Filter 


Exposure Time 




(J2000) 


(J2000) 




(s) 


NGC 541 


01:25:44.3 


-01:22:46 


F450W 


2x1200 










1x1000 








F814W 


2x 600 


NGC 2832 


09:19:46.9 


+33:44:59 


F450W 


4x1300 








F814W 


2x1300 


NGC 4839 


12:57:24.2 


+27:29:54 


F450W 


6x1300 








F814W 


2x1300 


NGC 7768 


23:50:58.6 


+27:08:51 


F450W 


6x1300 










1x1200 








F814W 


2 X 1300 










1x1200 




rAisLi: 2. Pro 


I'EHTIES OF 


lIo^T Cl.I 




Galaxy Cluster (cz)^ 






Lx" 


Richness Class'' BM Class'' 


(km s^-*-) 


(km s^-*-) 


(kcV) 


(lO*-* erg 


s-l) 



N541 


A194 


5388±61 


480+« 


2.63+0.15 


0.22+0.04 





II 


N2832 


A779 


6887+110 


503i^°° 




0.28+0.07 





I-II 


N4839 


A1656 


7339 


329 


4.5 








N7768 


A2666 


8086+105 


47g+95 




0.05 





I 



^ {cz) and cr from Zabludoflf, Geller & Huchra (1990), except for NGC 4839 (CoUess & Dunn 1996). 
X-ray temperatures and luminosities for NGC 541 from Wu, Xue &: Fang (1999). The X-ray temperature of NGC 4839 is taken from 
Neumann et al. (2003). 
" Richness and Bautz-Morgan classifications from Leir & van den Bergh (1977). 



Table 3. Surface Brightness Profiles 





Ms 




e 


PA 




(arcsec) 


(mag arcsec"^) 


(mag arcsec"^) 




(deg) 








NGC 541 








0.16 


15.507 


17.998 


0.07 


120.6 


-0.007 


0.28 


15.772 


18.292 


0.10 


128.6 


-0.004 


0.42 


16.012 


18.442 


0.07 


166.3 


0.018 


0.58 


16.175 


18.546 


0.01 


29.1 


-0.008 


0.76 


16.251 


18.682 


0.03 


72.8 


-0.005 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
Semi-major axis distance 
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Table 4. GC candidates photometry for NGC 541 



ID 


r 


A a 


AS 


B 


(B-I) 




(arcsec) 


(arcsec) 


(arcsec) 


(mag) 


(mag) 


1 


4.047 


-4.012 


-0.540 


24.200 


2.005 


2 


4.523 


0.087 


4.522 


25.102 


1.711 


3 


5.292 


-5.029 


-1.650 


24.545 


1.675 


4 


5.489 


-5.332 


-1.308 


25.594 


1.956 


5 


5.510 


-5.369 


1.243 


25.321 


1.861 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



Table 5. GC candidates photometry for NGC 2832 



ID 


r 


Aa 


AS 


B 


{B-I) 




(arcsec) 


(arcsec) 


(arcsec) 


(mag) 


(mag) 


1 


3.756 


-1.916 


3.402 


25.260 


2.113 


2 


4.117 


2.033 


-3.753 


25.842 


2.273 


3 


4.865 


4.644 


-2.960 


24.9.36 


1.813 


4 


4.941 


-2.327 


-4.546 


25.486 


1.653 


5 


5.286 


1.442 


5.149 


24.968 


2.113 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



Table 6. GC candidates photometry for NGC 4839 



ID 


r 


Aa 


AS 


B 


(B-I) 




(aix-.c-c) 




(arrs<-r) 


(mag) 


(mag) 


1 


3.969 


3.928 


1.900 


25.460 


2.194 


2 


5.818 


-0.369 


-5.809 


25.387 


1.565 


3 


6.070 


6.118 


2.718 


26.194 


1.883 


4 


6.182 


5.539 


3.752 


25.980 


2.121 


5 


6.277 


-6.548 


2.381 


25.137 


1.402 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



Table 7. GC candidates photometry for NGC 7768 



ID 


r 


Aa 


AS 


B 


(B-I) 




(arcsec) 


(arcsec) 


(arcsec) 


(mag) 


(mag) 


1 


3.593 


-2.871 


2.526 


25.824 


1.670 


2 


4.855 


1.898 


-4.552 


26.172 


2.120 


3 


5.717 


5.231 


3.319 


26.242 


2.038 


4 


8.337 


9.357 


-0.428 


25.635 


1.853 


5 


10.016 


1.635 


9.910 


26.362 


1.361 



Note. — 



The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
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Table 8. K-corrections 



lemplate 


[hc/ti\ 


cz 


z 


K(B) 


K(I) 




(dex) 


(km s j 




(mag) 


(mag) 




~ 




01 80Q + n 00002 


0.085 


0.020 








0231 8 + 00002 


0.111 


0.025 






7Qf;9_|_i 9 




0.118 


0.027 






Q1 QO-UI fi 




0.130 


0.031 


G2 


—0.4 




n 01 8nQ + n 00002 


0.072 


0.013 








n 02^^!^ 4- 00002 
u.u.^o±o 31 u.uuuu^ 


0.092 


0.017 








024^(S + 00004 


0.097 


0.018 








027*^2 -k OOOO'^ 


0.108 


0.021 


G3 


— 1.0 




01K0Q 4- 00002 
u.uiouy _L- u.uuuu^ 


0.052 


0.006 








02^^! 8 + 00002 


0.067 


0.007 






7362±12 


0.02456 ± 0.00004 


0.071 


0.007 






8192±16 


0.02732 ± 0.00005 


0.078 


0.008 


G4 


-1.5 


5422± 6 


0.01809 ± 0.00002 


0.0'M 


0.004 






6948± 6 


0.02318 ± 0.00002 


0.056 


0.006 






7362±12 


0.02456 ± 0.00004 


0.059 


0.006 






8192±16 


0.02732 ± 0.00005 


0.066 


0.007 


G5 


-1.9 


5422± 6 


0.01809 ± 0.00002 


0.038 


0.003 






6948± 6 


0.02318 ± 0.00002 


0.048 


0.004 






7362±12 


0.02456 ± 0.00004 


0.050 


0.004 






8192±16 


0.02732 ± 0.00005 


0.057 


0.005 



Table 9. K-correction Coefficients for (B — I)o 



Galaxy C'k Dk 

(mag) 



NGC 541 
NGC 2832 
NGC 4839 
NGC 7768 



0.96 0.032 

0.94 0.049 

0.94 0.056 

0.93 0.059 



Table 10. Statistics of the Globular Cluster Color Distributions 



Galaxy 






A7<= 






(mag) 


(mag) 







NGC 541 1.8ll°-»? 0.28t»-°} 0.2A+_°at 1-021^5 

NGC 2832 i.79t;;f2 o.28i;'/;;i o.i5i;!-f,i i.02i;!-;;^ 

NGC 4839 1.71«-°? 0.28i«:»i 0.21«-'S 1.0lt'^-°I 

NGC 7768 1.72+°;[;^ 0.28t^;'(!^ 0.47lj;;g 0.95t^;J| 



^Cbi = biwcight location estimator (Beers, Flynn & Gebhardt 1990). 
^Sbi = biweight scale estimator (Beers, Flynn & Gebhardt 1990). 
'^AI = assymetry index (Finch 1977). 
•^TI = tail index (HoagUn, Mosteller & Tukey 1983). 



Table 11. Observed and Derived Photometric Properties of cD Galax- 
ies 



Galaxy E(B - 1) {B - If^f - I)^;^ Vtot cz My.tot L/L*'' 





(mag) 


(mag) 


(mag) 


(mag) 


(km s ^) 


(mag) 




NGC 541 


0.080 


2.14 


1.81 


12.08±0.15 


5422± 6 


-22.44±0.28 


3.9 


NGC 2832 


0.039 


2.15 


1.79 


11.87±0.13 


6948± 6 


-23.10±0.26 


7.2 


NGC 4839 


0.023 


2.17 


1.71 


12.06±0.10 


7362±12 


-23.02±0.25 


6.7 


NGC 7768 


0.087 


2.22 


1.72 


12.30±0.13 


8191±16 


-23.10±0.26 


7.2 



aPor Ma = -21.79 + 51og/t65 (Mobasher et al. 2003), Hq = 72 km s"! Mpc"! (Preedman et al. (2001), and (V - R) = 0.61 Pukugita, 
Shimasaku & Ichikawa (1995). 
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Table 12. GLonri.AR Ci,rsTEB Specific Frech'ekcies 

Galaxy M™<=t ATg-g* S'^"^ 
(mag) 

NGC 541 -21.40 ±0.2 1620±450 4.4 ± 1.5 

NGC 2832 -22.35 ± 0.2 3200±920 3.7 ± 1.3 

NGC 4839 -21.84 ±0.2 3060±850 5.6 ± 1.9 

NGC 7768 -22,:«±n,2 1850±570 2,2 ±0.8 



